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ABSTRACT 

The luminosities and effective temperatures, as well as the whole bolometric 
lightcurves of nonlinear convective RR Lyrae models with O.OOOl^Z^O.006 are trans- 
formed into the SDSS photometric system. The obtained ugriz lightcurves, mean mag- 
nitudes and colors, pulsation amplitudes and color-color loops are shown and analytical 
relations connecting pulsational to intrinsic stellar parameters, similarly to the ones 
currently used in the Johnson-Cousins filters, are derived. Finally the behaviour in 
the color-color planes is compared with available observations in the literature and 
possible systematic uncertainties affecting this comparison are discussed. 
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Among Population II radially pulsating stars, RR Lyrae 
play a relevant role both as standard candles and as stel- 
lar property tracers. They are bright stars, easily identi- 
fied thanks to their characteristic variability (light curves, 
periods, and colors) and their luminosity span a narrow 
range. They usually show low {Z ~0.0001) to intermedi- 
ate [Z ~ 0.006) metallicity, but they can reach solar abun- 
dances in the solar neighbourhood or in the Galactic bulge. 
The calibration of their absolute magnitude Mv{RR) in 
terms of the measured iron-to-hydrogen content [Fe/H] al- 
lows us to use these variables to infer the distances of Galac- 
tic Globular Clusters (GGC) and nearby galaxies, as well as 
to calibrate secondary d istance indicators such as the GC 
luminosity function (see iDi Criscienzo et alJ l2006l and ref- 
erences therein). In this context several observational and 
theoretical efforts have been made in the last years to pro- 
vide accurate evaluations of the Mv(RR)-[Fe/H] relation, 
usu ally approximated as Mv(RR) = a + f5[Fe/H] (see 
e.g. ICaputo et allbOOl ICacciari et alll2003l and references 
therein), and to derive relevant implications for the Pop. II 
distance scale and the age of globular clusters. The above 
characteristics make RR Lyrae superb probes of the old 
stellar populations and studies of their pulsation properties 
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have provided much of our present knowledge of the struc- 
ture, kinematics, and the metal abundance distribution of 
the halo. Moreover, they have been adopted in the recent 
literature as fundamental targets of surveys devoted to the 
identification of specific populations and galactic substruc- 
tures (see e.g. lWu et alJ 2005. Brown ot al. 2004, Vivas ct all 
I2OO4I and references therein). In pa rticular, the Sloan Dig- 
ital Survey (SDSS) for RR Lyrae dlvezic et al.ll2000tl and 
the QUEST RR Lyrae survey (Vivas et al. 2001, 2004, 
2006) detected the tidal stream from the Sagittarius dSph 
galaxy and other density enhancements in the halo that 
may be other tidal streams, supporting the idea that RR 
Lyrae surveys are crucial to trace the merger history of the 
Milky Way. In t he context of the VLT Survey Telescope 
(VST) GTO (see lAlcala et al."2006^ we have planned a sur- 
vey (STREGAQVST, see Marconi et al. 2006) devoted at 
the exploration of the southern part o f the Fornax stream 
jLvnden-Bellll982LlDinescu et alJl2004l) and the tidal inter- 
action of the involved satellite galaxies and globular clus- 
ters with the Milky Way halo. To this purpose RR Lyrae 
will be used as tracers of the oldest stellar populations 
through multi-epoch observations in the SDSS g and i filters. 
The comparison with theoretical period-luminosity-color, 
peri od-luminosity-amplitude and Wesenheit relations (see 
e.g. iMarconi et al.ll20o3 : | Di Criscienzo. Marconi fc Caputol 
I2OO4I) will provide information on the individual distances 
and in turn on the spatial distribution of the investigated 
stellar system. In this context, the addition of multi-epoch 
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Table 1. Model input parameters. 



z 


Y 


M/Mq 


log L/ log Lq 


0.0001 


0.24 


0.65 


1.61 


0.0001 


0.24 


0.70 


1.72 


0.0001 


0.24 


0.75 


1.61 


0.0001 


0.24 


0.75 


1.72 


0.0001 


0.24 


0.75 


1.81 


0.0001 


0.24 


0.80 


1.72 


0.0001 


0.24 


0.80 


1.81 


0.0001 


0.24 


0.80 


1.91 


0.0004 


0.24 


0.70 


1.61 


0.0004 


0.24 


0.70 


1.72 


0.0004 


0.24 


0.70 


1.81 


0.001 


0.24 


0.65 


1.51 


0.001 


0.24 


0.65 


1.61 


0.001 


0.24 


0.65 


1.72 


0.001 


0.24 


0.75 


1.61 


0.006 


0.26 


0.58 


1.55 


0.006 


0.26 


0.58 


1.65 


0.006 


0.26 


0.58 


1.75 



r and single epoch u exposures on selected fileds will enable 
us to derive further constraints on the intrinsic stellar pa- 
rameters of RR-Lyrae stars. However, in order to correctly 
compare model predictions with observations all the pre- 
dicted pulsation observables need to be transformed into 
the SDSS photometric filters. To this purpose in this paper 
we use updated model atmospheres to predict the pulsation 
properties of RR Lyrae stars of different metal contents in 
the u, g, r, i, z bands^ . 

The organization of the paper is the following: in Sec- 
tion 2 we present the adopted pulsation models, in Section 
3 we discuss the procedures adopted to transform the the- 
oretical scenario into the SDSS filters and in Section 4 we 
illustrate the pulsation observables in these bands. Finally 
in Section 5 we compare model predictions with SDSS RR 
Lyrae data available in the literature, discussing possible 
systematic errors affecting the comparison. The Conclusions 
close the paper. 



2 THE ADOPTED PULSATION MODELS 

During the last few years we have been computing an 
extensive and detailed set of nonlinear nonlocal time- 
dependent convective models for RR Lyrae stars, spanning 
a wide range of p hysical parameter s and chemical com- 
positions (see e. g. |Bono_et_alJ E oOSl. iMarconi et all l2003l 
iDi Criscienzo. MMCon^^^aputol^20M) . In this paper we 
concentrate on models with metallicity between Z = 0.0001 
and Z = 0.006, that is the typical range for Galactic 
Globular Cluster RR Lyrae. The adopted metallicity. He- 
lium content and stellar parameters are reported in Ta- 
ble 1, whereas the pulsation properties of these models in 
the Johnson-Cousins bands, a s obtained by adopting the 
static model atmospheres by ICastelli. Gratton fc Kurucj 
1^93), as well as an extensive comparison with observed 
RR Lyrae, are discussed in our previous papers (see e.g. 

^ These are the filters that will be mounted on the VST. 



Di Criscienzo. Marconi fc Caputd 12004 hereinafter D04, 
Bono et alJl200A iMarconi et alJl2003llBono et al.lll997ri . In 

particular D04 discuss in detail the dependence of pulsa- 
tion properties on the adopted mixing lenght parameter 
(q = I /Hp where I is the mixing lenght and Hp the pres- 
sure height scale) that enter the turbolent-convective model 
to close the n onlinear system of dynami cal and convective 
equations (see lBono fc Stellingwer f 1994Y From this analy- 
sis the authors conclude that the standard value {a = 1.5) 
adopted in all their previous investigations, well reproduces 
the behaviour of RR Lyrae in the blue regions of the insta- 
bility strip, whereas there are indications for an increasing a 
value (up to 2.0) as one moves toward the red edge. In this 
paper we concentrate on the models computed with a = 1.5 
and the stellar parameters reported in Table 1, but we will 
mention the effect, if any, of a possible alpha increase. 



3 TRANSFORMATION INTO THE SDSS 
PHOTOMETRIC SYSTEM 

In order to obtain the pulsation observables of the investi- 
gated RR Lyrae models in the SDSS bands we have trans- 
formed both the individual static luminosities and effective 
temperatures and the predicted bolometric light curves into 
the corresponding filters. In particular, as the adoption of 
linear transformations from one photometric system to an- 
other introduces uncertainties (limited precision and strong 
dependence on the color range), we directly build magni- 
tudes and colors for our RR-Lyrae models by convolving 
model atmosphere fluxes with SDSS transmission functions 
^ . In general, the calculation of the magnitudes i n a given 
photometric system (see e.g. iGirardi et aljl2'002h involves 
the integral equation: 

ms, = -2.5 log -^1— + mo (1) 

where Sx is the transmission function, f\ is the stellar 
flux (that corresponds to model atmospheres of known 
(Te//, [M/_ff], log (/)), f1 the zero point reference flux and 
mo the zero point reference magnitude. Then, the final ab- 
solute magnitudes are computed by the knowledge of the 
stellar radius. 

In the SDSS photometric system (an ABmag system), 
mo = and the zero poi nt reference spe ctrum is the absolute 
flux of Vega at 5480 A iFukugitalll99^ . 

As for model atmospheres, in this paper we 
adopt the homogeneous set of updated ATLAS9 
Kurucz model atmospheres and synthetic fluxes 
(new-ODF models)"^ computed with a new set of 
Opacity Distribution Functions iCastelli fc Kurric3 
(2003) • These calculations assume steady-state plane- 
parallel layer atmospheres, covering a metallicity range 
[M/H] : 0.5, 0.2, 0.0, -0.5, -1.0, -1.5, -2.0 and - 2.5, both 
for [a/Fe] = 0.0 and [ct/Fe] = 0.4, temperatures between 
3500 to 50000 K and log 3 from 0.0 to 5.0. The lKuruc3 

^ The SDSS transmission curves are available at url 
jhttp: / /www. sdss.org/dr3/instruments/imager/index. html 

Available at url http://kurucz.harvard.edu/grids.html or 
pittp: / / wwwuser . oat .ts. astro, it / castelli /grids . html| 
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l|l990(l model atmospheres had some recognized problems, 
in particular for effective temperatures lower than 4500 
K. One of the reasons was the lack of H2O in the line 
opacity calculations and the use of approximate line data 
for TiO and CN. The ATLAS9 grids of the new-ODF 
mod els were computed with updated solar abundances 
from iGrevesse fc Sauvall (^9^ and considering molecular 
line lists which include H2O molecular transitions and 
updated TiO and CN data. A comparison between broad- 
band synthetic colors for late type giants computed from 
the new-ODF models and from the PHOENIX/NextGen 
models has shown a remarkable agreement (Kucinskas et 
al. 2005, 2006) in spite of the PHOENIX/NextGen models 
are computed with more molecular species than in ATLAS9 
and by assuming spherical simmetry. This last hypothesis 
is very important when the extent of the atmosphere is 
comparable with the radius of the star. RR-Lyrae stars 
have generally temperatures higher than 4000 K, however, 
as a sanity check, for selected models we also computed 
magnitudes and colors by convolving the light curves with 
the PHOENIX/NextGen synthetic spectra * (see Sect. 5). 



4 THE PULSATION OBSERVABLES IN THE 
SDSS PHOTOMETRIC SYSTEM 

The application of the procedure discussed in the previous 
Section allows us to predict the pulsation observables and 
their behaviour as a function of the model input parameters 
in the SDSS filters. In the following we explore in detail the 
morphological features of light curves and the behaviour of 
mean magnitudes and colors, the topology of the instability 
strip and color-color loops, as well as the main relations 
between the periods (and the amplitudes) of pulsation and 
the intrinsic stellar parameters. 

4.1 The light curves 

All the computed bolometric light curves have been trans- 
formed into the SDSS bands ^. Figs. [TO show the trans- 
formed light curves for both fundamental (F) and first 
overtone (FO) pulsators at selected luminosity levels and 
Z = O.OOOi, 0.001, 0.006. 

As already found for the Johnson- Cousins bands, the 
pulsation amplitudes vary with the wavelenght, increasing 
from utog and decreasing from g to z. Moreover, for a fixed 
filter, the amplitudes of fundamental pulsators decrease from 
the blue to the red edge, whereas the first overtone ones in- 
crease moving from the blue edge to the middle of the in- 
stability strip and decrease as the red edge is approached. 
The agreement between these predicted trends and the ob- 
served behaviour of both fundamental and first overtone pul- 
sators has been extens ively discussed in previous pap e rs (see 
e.g. the discussi o n in iBrocato. Castellani fc Ripepil ll 19961) 
and iBono et alJ lll997ll '). We also notice that the theoret- 
ical light curves in the Johnson-Cousins bands have been 
successfully compared with available data in the literature 



(see lBono. Castellani fc Marconil (l2000h . ?). Unfortunately a 
similar comparison cannot be performed for the SDSS filters 
because of the lack of well-sampled RR Lyrae light curves 
in these bands. 



4.2 The mean magnitudes and colors 

From the obtained u, g, r, i, z light curves we can de- 
rive mean magnitudes and colors, following different aver- 
aging procedures to produce either magnitude-averaged or 
intensity-averaged values. As already extensively discussed 
for the Johnson- Cousins bands (see D04) the various types 
of mean values differ from the static ones the stars would 
have were they not pulsating. In Fig. l4l5l we show (for F and 
FO-models respectively) the difference between magnitude- 
averaged and intensity-averaged mean magnitudes for the u, 
g, r, i, z filters (top panels), as well as the differences be- 
tween the two kinds of mean values and the corresponding 
static quantities (middle and bottom panels). 

We notice that similarly to what happens for the 
Johnson- Cousins filters, the difference between magnitude- 
averaged and intensity-averaged magnitudes increases as the 
corresponding pulsation amplitude increases, reaching the 
highest values in the u and g bands. Moreover we confirm 
that intensity-averaged magnitudes better reproduce the be- 
haviour of static values than the magnitude-averaged ones. 
Indeed the difference between the latter mean values and 
the static ones can reach 0.2 mag in the g band. For this 
reason in the following we will adopt the intensity-averaged 
mean magnitudes in the u, g, r, i, z filters. These quantities 
are reported in Table 2 and 3 for the whole fundamental and 
first overtone model sets respectively®. 

4.3 The instability strip boundaries 

In Fig.|S|we show the predicted instability strip boundaries 
for both fundamental (solid lines) and first overtone (dashed 
lines) models in the intensity-averaged g versus g-r plane, for 
the metal abundances labelled in the different panels. From 
left to right the different lines indicate the first overtone 
blue edge (FOBE), the fundamental blue edge (FBE), the 
first overtone red edge (FORE) and the fundamental red 
edge (FRE). At each magnitude level blue (red) boundaries 
correspond to the effective temperatures of the first (last) 
pulsating model in the selected mode. Concerning the be- 
haviour of the boundaries in the magnitude-period diagram, 
a linear regression to all the models quoted in Table 1 yields, 
by taking into account the first (for the FOBE) and the last 
(for the FRE) pulsating models, the mass-dependent ana- 
lytical relations reported in Table 4. The zero point of the 
FOBE relation vary by ~0.08 mag in g and by ~0.06 mag 
in i if the assumed a value increases from 1.5 to 2.0. The 
corresponding zero-point variation for the FRE is ~0.3 mag 
in g and ~0.01 mag in i (see also the discussion in D04). 

Only for the u filter we also find a non negligible de- 
pendence on the metal content and in all cases the standard 
deviations are of the order of few hundredths of magnitude. 



* Available at the web site ftp: / /ftp.hs.uni-hamburg.de/pub/outgoing/phoenix /GAIA/' 

^ The full set of transformed ligiit curves is available upon request ° Similar tables for magnitude-averaged values are available upon 

to the authors. request to the authors. 



4 M. Marconi, M. Cignoni, M. Di Criscienzo, V. Ripepi, F. Castelli, I. Musella, A. Ruoppo 




Predicted properties of RR Lyrae stars in the SDSS photometric system 5 




10 10 10 10 1 

Phase 

Figure 2. Predicted SDSSS light curves for selected F (top) anf FO (bottom) models with Z=0.001, M = O.65M0,log L/ log = 1.61 
and varying eflfective temperature. 
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10 10 10 10 1 

Phase 

Figure 3. Predicted SDSSS light curves for selected F (top) anf FO (bottom) models with Z=0.006, M = O.58M0,logL/logL0 = 1.65 
and varying effective temperature. 
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Table 2. Intensity-averaged mean magnitudes for the full set of fundamental models (the full table is available in the electronic form). 



Z M/Mq log L/ log Lq Te (K) P (d) <u> <g> <r> <i> <z> 



0.0001 0.65 1.61 6900 0.4059 1.8978 0.8378 0.7813 0.8054 0.8419 

0.0001 0.65 1.61 6800 0.4260 1.9090 0.8445 0.7675 0.7808 0.8118 



Table 3. Intensity-averaged mean magnitudes for the full set of first overtone models (the full table is available in the electronic form). 



z 


M/Mq 


log L/ log Lq 


Te (K) 


P (d) 


<u> 


<9> 


< r > 


< i > 


< z > 


0.0001 


0.65 


1.61 


7300 


0.2527 


1.8838 


0.7553 


0.7572 


0.8255 


0.8918 


0.0001 


0.65 


1.61 


7200 


0.2634 


1.8843 


0.7737 


0.7664 


0.8242 


0.8826 
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Figure 4. Comparison between static and mean magnitudes for 
F-models. 



Figure 5. The same as the previous figure but for FO-models. 



Table 4. Analytical relations for the FOBE and FRE in the form 



Mi = a + blogP + cl 
also reported. 



M 



+ dlogZ. The dispersion a (mag) is 





boundary 


a 


b 


c 


d 


tj 


u 


FOBE 


0.40 


-2.32 


-1.95 


0.086 


0.03 




FRE 


1.64 


-1.79 


-2.43 


0.12 


0.03 


g 


FOBE 


-1.04 


-2.43 


-2.14 




0.04 




FRE 


0.18 


-1.99 


-2.34 




0.02 


r 


FOBE 


-1.12 


2.63 


-1.85 




0.03 




FRE 


-0.07 


-2.23 


-1.92 




0.02 


i 


FOBE 


-1.11 


-2.72 


-1.86 




0.03 




FRE 


-0.18 


-2.34 


-1.96 




0.01 


z 


FOBE 


-1.09 


-2.78 


-1.91 




0.02 




FRE 


-0.22 


-2.40 


-2.05 




0.01 



4.4 The color-color loops 

The transformed multifilter light curves can be reported in 
diflterent color-color planes. In Figs. l7l9l we show the theoret- 
ical loops in the g-r vs u-g,r-i vs g-r and i-z vs r-i diagrams 
for selected models (see Table 5) at each adopted metallicity. 

We notice that the metallicity effect is more evident in 
Fig. [7| due to the significant sensitivity of the u band on 
metal abundance. On the other hand, in the r — i vs g — r 
plane the loops are very narrow and the effect of metallic- 
ity is smaller, so that, once the metallicity is known, the 
comparison between theory and observations in this plane 
could be used to evaluate color excesses. As for the i — z 
vs r — i diagram, the predicted loops are very close to each 
other and the metallicity dependence is much less evident 
than for the other color combinations. A linear regression 
through the intensity-averaged mean magnitudes reported 



8 M. Marconi, M. Cignoni, M. Di Criscienzo, V. Ripepi, F. Castelli, I. Musella, A. Ruoppo 




Figure 6. Predicted instability strip boundaries for both F (solid 
lines) and FO (dashed lines) models in the intensity-averaged 
g versus g-r plane. The adopted metal abundances and stellar 
masses are labelld. 
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Figure 7. Theoretical loops in the g-r vs u-g diagram for selected 
models (see Table 5) at each adopted metallicity. 



Table 5. Physical parameters of models plotted in the color-color 
diagrams. 
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Y 


M/Mq 


log L/ log Lq 




mode 


0.0001 


0.24 


0.80 


1.81 


7100 


FO 


0.0001 


0.24 


0.80 


1.81 


6800 


FO 


0.0001 


0.24 


0.80 


1.81 


6700 


F 


0.0001 


0.24 


0.80 


1.81 


6300 


F 


0.0001 


0.24 


0.80 


1.81 


5900 


F 


0.0004 


0.24 


0.70 


1.72 


7000 


FO 


0.0004 


0.24 


0.70 


1.72 


6700 


FO 


0.0004 


0.24 


0.70 


1.72 


6700 
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0.0004 


0.24 


0.70 


1.72 


6400 
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0.0004 


0.24 


0.70 


1.72 
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0.001 


0.24 


0.65 


1.61 


7200 


FO 
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1.61 
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0.65 


1.61 


6700 
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0.24 


0.65 


1.61 
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0.001 


0.24 


0.65 


1.61 


6100 


F 


0.006 


0.26 


0.58 


1.65 


7100 


FO 


0.006 


0.26 


0.58 


1.65 


6900 


FO 


0.006 


0.26 


0.58 


1.65 


6800 


FO 


0.006 


0.26 


0.58 


1.65 


7000 


F 


0.006 


0.26 


0.58 


1.65 


6400 


F 


0.006 


0.26 


0.58 


1.65 


6000 


F 



in Tables 2 and 3 provides the following metal-dependent 
analytical color-color relations: 

?t-5 = 1.41 -0.12(3 -r) + 0.088 log Z (a = 0.03) (2) 



g-r = 0.234 + 2.11(r - i) + 0.035 log Z {a ^ 0.006) 

r - i = 0.085 -I- 2.04(j ~ z) + 0.0097 log Z {a ^ 0.004) 

holding for fundamental pulsators, and: 

u~-g = 1.35 - 0.19(g - r) + 0.065 log Z {a = 0.02) 

g~r = 0.217 + 1.95(r - i) + 0.0258 log Z (a ^ 0.006) 
r~i = 0.082 + 2.00(j - z) 

for the first overtone ones. 



(3) 



0.0062 log ^ ((7 = 0.003) 
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Figure 8. The same as in Fig.lTlbut in the r — ivsg — r diagram; 
in this plane the loops move upward as the metallicity decreases. 
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Figure 9. The same as in Fig.|7|but in the i — z vs r — i diagram 
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4.5 The multifilter Period-Magnitude-Color 
relations 

The obvious outcome of the pulsation relation (see van Al- 
bada & Baker 1971, Di Criscienzo et al. 2004) connecting the 
period to the mass, the luminosity and the effective temper- 
ature into the observational plane is the Period-Magnitude- 
Color (PMC) relation, where the pulsation period for each 
given mass is correlated with the pulsator absolute magni- 
tude and color. A linear regression through the intensity- 
averaged values reported in Tables 2 and 3 provides a PMC 
relation for each selected couple of bands and chemical com- 
position. 

In particular we find : 

u = 6.43 - 0.89 log P - 3.78(« - g) - 2.841og Af -f 0.36 logZ 
a = 0.14 (4) 
g = -1.05 - 2.871ogP-f 3.35(g - r) - 1.87 log J\/ - 0.06 logZ 
o- = 0.02 (5) 
g = -0.83 - 2.871ogP-f 2.27(g - i) - 1.94 logM - 0.024 logZ 
a = 0.02 (6) 
g = -0.67- 2.91 log P + 1.99(3 - 2) - 1.90 log M 
cr = 0.07 (7) 

for fundamental models and 

u = 5.77- 1.48 log P - 3.87(« - g) - 2.72 log M -f 0.25 logZ 
o- = 0.07 (8) 
g = -1.53 - 3.111ogP + 3.57(g - r) - 1.63 log M - 0.042 log Z 
cr = 0.01 (9) 
g = -1.279 - 3.071ogP-|-2.334(g - i) - 1.72 log M - 0.017 log Z 
a = 0.009 (10) 
g = -1.151 - 3m9logP + 2.02(g - z) - 1.77 log M - 0.0045 log Z 
o- = 0.007 (11) 

for the first overtone ones. As discussed in D04 for the Johnson- 
Cousins bands, these relations allow us to derive an estimate of 
the stellar mass for RR Lyrae of known distance and color, while 
for cluster pulsators sharing the same distance and reddening, 
they provide direct estimates of the mass spread. On the other 
hand, if the mass and colors are known, the same relations can 
be used to infer individual and/or mean distance moduli of RR 
Lyrae stars in a given globular cluster or galaxy. 



4.6 The pulsation amplitudes 

The pulsation amplitudes of the fundamental and first overtone 
multiwavelenght lightcurves are reported in Tables 7 and 8. Their 
behaviour as a function of the pulsation period is shown in Fig. 
1101 for selected model sequences at Z = 0.0001 (left panels) and 
Z = 0.001 (right panels). We notice that the fundamental pul- 
sators follow a linear behaviour, thus allowing the derivation of 
a mass dependent linear relation between the period, the pulsa- 
tion amplitude and the absolute magnitude. The coefficients of 
these relations for the whole fundamental model set and the var- 
ious photometric bands are reported in Table 6. We notice that 
the coefiicients of these relations are expected to depend on the 
adopted a parameter (see also D04). In particular in the g filter 
the zero point and the amplitude coefficient should vary by ~ 
0.14 and ~ 0.38 respectively, as a increases from 1.5 to 2.0. 

As for first overtone pulsators we notice the characteristic 
bell shape also shown in the bolomctric and Johnson-Cousins 
band Bailey diagram (see e.g. Bono ct al. 1997). In Figs. ITT1 and 
I12l we show the amplitude ratios between the Johnson V and the 



Table 6. Analytical coefficients of the PL A relations in the form 



logP : 



- bA, +c< Mi> +dlog -I- elog Z. 





a 


b 


c 


d 


e 


cr 


u 


0.72 


-0.186 


-0.41 


-0.73 


0.038 


0.02 


g 


0.14 


-0.184 


-0.371 


-0.54 


0.009 


0.02 


r 


-0.03 


-0.174 


-0.372 


-0.60 


-0.015 


0.01 


i 


-0.07 


-0.171 


-0.370 


-0.64 


0.069 


0.01 


z 


-0.074 


-0.164 


-0.368 


-0.66 


0.074 


0.009 



SDSS g bands (left top panel of each figure) and between the u, 
r, i, z and the g bands, for fundamental and first overtone models 
respectively. We notice that the g band amplitude is systemati- 
cally higher than the V amplitude, independently of the period 
and the adopted metal abundance (see labels). At the same time 
the amplitudes in the r,i, z filters scale with a constant mean ratio 
(ranging from about 0.7 to 0.5 from r to z) with the g band am- 
plitude, thus suggesting that only few points along the lightcurves 
in these filters will be required, if the g curve is accurately sam- 
pled. As for the u band, the scatter is larger due to the significant 
dependence on the adopted metallicity. This occurrence is more 
evident for fundamental models, which are characterized by more 
asymmetric light curves and higher pulsation amplitudes. 

We note that, once known the metal content, the mass term 
contained in all the analytical relations reported in this Section 
can be replaced with the value predicted by evolutionary hori- 
zontal branch computations at the selected me tallicity (see e.g. 
columns 5 and 6 of Table 1 in iBono et alJl2003^ . 



5 COMPARISON WITH THE OBSERVATIONS 

In this section the predicted color-color loops in the various filter 
combinations are compared with observed RR Lyrae samples. In 
par ticular we test ou r predictions with the QUEST survey (see 
e.g. *^^va^^al.l[2 0pi ]) and the observations in the Draco galaxy 
iBo nanos et al.l l2004lh. The first one, a large census of field RR 
Lyrae, should map the history of the outer halo of our Galaxy, 
the second should represent an example of RR Lyrae in a difl'erent 
enviroment. 

Both the samples are positionally matched against the SDSS- 
DR4 iAdelman-McCarthv et"ai]|2006^ catalogue using a search 
ra dius of 0.1 arcsec and subsequ ently d ereddened (using IR maps, 
see lSchlegei . Finkbeiner & Dayi sll993) . Figs.ll3landll4lshow the 
distributions for QUEST data, as observed in the g-r vs u-g and 
r-i vs g-r plane respectively. Superimposed, we show the theoret- 
ical loops for the labeled metallicities instead of average values 
because SDSS observations consist of one or, at most, two phase 
points. The QUEST data which are expec ted to trace metal poor 
(Z < 0.001) and distant halo stars (see IVivas et al.ll20oIh are 
not matched by the corresponding model loops, with the discrep- 
ancy larger than the mean uncertainty resulting from photomet- 



ric (cr(u — g) < 0.03 ma g) and reddening errors (less than 0.01 
mag in colors, see also llvezic et all l200fi) . In Fig. 1151 we show 
the comparison with DRACO RR Lyrae in the g-r vs u-g plane. 
In this case we have a large spread in u-g due in part to the 
significant photometric uncertainties (cr(u — g) ~ 0.1 mag and 
""(9 — r) 0.02) affecting the data. In this case, the compari- 
son with the theoretical loops does not allow us to discriminate a 
metallicity effect. However, the mean metallicity of this galaxy i s 
generally considered poorer than Z=0.001 (see e.g. lMateolll99d) . 
thus it is noteworty that a consistent fraction of the Draco RR 
Lyare is redder than our Z=0.006 models with an overdensity at 
u-g > 1.25 (difficult to explain with the photometric error alone). 
The situation is similar in the r-i vs g-r plane (see Fig. 
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Table 7. Pulsation amplitudes for the full set of fundamental models (the full table is available in the electronic form). 



z 


M/Mq 


log L/ log Lq 


Te (K) 


P (d) 




A, 






A, 


0.0001 
0.0001 


0.65 
0.65 


1.61 
1.61 


6900 
6800 


0.4059 
0.4260 


1.7877 
1.6424 


1.8675 
1.7906 


1.3689 
1.2955 


1.0900 
1.0261 


0.9669 
0.9187 




Figure 10. Period-amplitude diagram for FO (dashed line) and F (solid line) models with M = 080AfQand Z=0.0001(left panel) and 
M = O.65M0 Z=.001(right panel) and for the three labelled values of logL/logL0. In both cases star symbols represent models with 
M = 0.75Mq. 



ITell . Possible sources for the discrepancies between theory and ob- 
servation may result from observational uncertainties, including 
reddening and contamination effects, or theoretical biases, as for 
example the adopted model atmospheres and chemical mixture. 

As for the reddening, a critical point is represented by the ex- 
tinction law adopted to relate the measured E(B-V) to the extinc- 
tions in the SDSS filters (s ee e.g. [^ira^i et^.l l2004l . Assuming 
the coefficents tabulated bv lGirardi et al . U2004h . an underestima- 
tion of the E(B-V) color excesses of the order of 0.02 mag would 
require a redshift in the theoretical u-g and g-r of about 0.03 and 
0.02 respectively, in the direction of reducing the discrepancy be- 
tween data and metal-poor predictions in Figs. ll3l and ll4l (see the 
arrow in these plots). 

Concerning model atmosp heres, as discussed in sect ion 3, we 
have adopted ATLAS models iCastelli fc KurucJl2003tl without 
a enhancement. In order to explore the effect of possible a en- 
hancement for the lowest metallicities, as empirically sugjgeste d 
by various authors in the literature (see e... l G;atton et alM . 



we have transformed again the theoretical bolometric light curves 
by using ATLAS model atmospheres with [a/Fe] =0.4 and re- 
ducing the adopted \ Fe/H] in order to obtain the sam e global 
model metallicity (see ISalaris. Chieffi &: Stranierolll993l) . As an 
example, in Fig. ll7l we show the variation of the theoretical loops 
for models with Z=0.0004. It appears that increasing [a/Fe] from 
to 0.4 produces a larger discrepancy between metal-poor models 
and the data. On the other hand, if in the conversion from Z to 
[Fe/H] wc replac e the assumed solar metallic i ty wit h the more 
updated value bv lAsplund. Grevesse fc Sauvall i2005l) . loops of a 
given metallicity move redward by about 0.05 mag in g-r vs u- 
g plane (see e.g. figure [TSl for Z=0.001), so that the recovered 
global metallicity for halo RR Lyare is closer to typical values in 
the literature. 

To take into account the effect of the adopted set of model 
atmospheres, Fig. ll9l Dlots the same comparison as in Fig. ll3l but 
using the PHOENIX atmospheres for giant stars I Kucinskas ct al] 
l2006r . which have the main advantage of assuming the spherical 
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Figure 11. Amplitude ratios between the Johnson V and the SDSS g bands (left top panel) and between the u, r, i, z and the g bands 
for fundamental models. 

Table 8. Pulsation amplitudes for the full set of first overtone models (the full table is available in the electronic form). 



z 


M/Mq 


log L/ log Lq 


Te (K) 


P (d) 




A, 






A, 


0.0001 
0.0001 


0.65 
0.65 


1.61 
1.61 


7300 
7200 


0.2527 
0.2634 


0.3939 
0.9037 


0.4978 
1.0765 


0.3522 
0.7677 


0.2727 
0.5969 


0.2427 
0.5296 



geometry instead of the classical plane-parallel structure. From 
the u-g vs g-r plot, it is evident that the PHOENIX atmospheres 
affect the u-g color with a blueward shift (of about 0.03 mag) 
for models with g-r between 0.1 and 0.3 mag, corresponding to 
the bulk of the observed RR Lyrae. Moreover, the PHOENIX 
atmospheres lead to redder u-g colors (by about 0.05 mag) for 
stars with g-r lower than 0.1 mag, where the RR Lyare density is 
lower. However, both these effects produce minor changes in the 
behaviour observed in Figs. [T^ and ll4l 



are found for r,i, z with respect to the g band, suggesting that 
only few points along the lightcurves in these filters will be re- 
quired, if the g curve is accurately sampled. The theoretical g-r 
vs u-g and r-i vs g-r are compared with available data in the lit- 
erature. In particular, for field RR Lyrae stars from the QUEST 
survey the comparison seems to suggest a higher mean metallic- 
ity than usually assumed. In order to understand this occurrence 
several sources of systematic errors have been discussed. 



6 CONCLUSIONS 

We have transformed the predicted scenario for RR Lyrae stars 
with Z in the range 0.0001-0.006 into the SDSS photometric sys- 
tem, providing theoretical tools for the interpretation of modern 
large RR Lyrae surveys in these bands. Mean magnitudes and 
colors and the pulsation amplitudes are used to derive multiband 
analytical relations that can be used to constrain both the dis- 
tances and the stellar masses. Rather constant amplitude ratios 
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Figure 12. The same as the previous figure but for first overtone models 
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Figure 13. Comparison between the theoretical loops shown in 
Fig. Q and the QUEST RR Lyrae data (see text for details). 
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Figure 14. Comparison between the theoretical loops shown in 
Fig. Eland the QUEST RR Lyrae data (see text for details). 



to the "Progetto 6: Calcolo Evoluto e sue applicazioni (RSV6) 
Cluster Cll/B". 
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Figure 19. The QUEST RR Lyrae data compared with se- 
lectd models transformed both with ATLAS (solid lines) and 
PHOENIX (dotted lines) model atmospheres. 
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